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Abstract 

We study theoretically the lattice instability in the spin gap systems un- 
der magnetic field. With the magnetic field larger than a critical value h&, 
the spin gap is collapsed and the magnetization arises. We found that the 
lattice distortion occurs in the spindadder at an incommensurate wavevector 
corresponding to the magnetization, while it does not occur in the Haldane 
system. At low temperatures the magnetization curve shows a first order 
phase transition with this lattice distortion. 
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The one-dimensional (ID) quantum spin systems with the gap in the excitation spec- 
trum have been attracting intensive interest both theoretically and experimentally This is 
due to the underlying large quantum fluctuations which destroy the magnetic long range 
order (LRO) Another aspect of these spin gap systems is that they are free from the 
instabilities, e.g., spin-Peierls instability [0, which comes from the diverging generalized sus- 
ceptibility as the temperature T — >• 0. One of the spin gap systems of the current interest is 
the antiferromagnetically coupled two- leg spin ladder system §-0. Under a magnetic field 
at T = the magnetization remains zero until h reaches a critical value h c \ because the 
ground state is the spin singlet. At h — h c \ the singlet and lowest triplet states are degener- 
ate, and the magnetization begins to appear. Between h c % and h c2 , where the magnetization 
increases, the system remains gapless and is described as a Tomonaga-Luttinger liquid with 
the continuously changing exponent. For h > h c2 the system is perfectly spin polarized. 
Recently the two-leg ladder spin system Cu2(l,4-Diazacycloheptane)2Ci4 has been studied 
under magnetic field The magnetization process has been studied also theoreti- 

cally P^[~P]. In addition, an anomaly of the specific heat is observed at T = IK when 



h = 9.0T > h c \ = 7.2T |Q. This strongly suggests the phase transition to the ordered state, 
which is assumed to be AF long range ordering |J , but is not yet confirmed. 

Motivated by these experimental works, we study in this letter the alternative instability, 
i.e., lattice instability, of the two- leg spin ladder system in its gapless region h cl < h < h c2 . 
We found that when coupled to the lattice distortion, which should be always there, the gap- 
less state is unstable to the (generalized) spin-Peierls state accompanied with the modulation 
of the inter-chain exchange interaction with the incommensurate wavenumber proportional 
to the magnetization. Another consequence of this instability is that the magnetization 
curve at low temperature shows a first order phase transition, and the schematic phase dia- 
gram is given in Fig. 1. This is in sharp contrast to the Haldane system, where there is no 
such an instability. 

We start with the bosonized Hamiltonian for the two coupled spin 1/2 Heisenberg model 
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where tpj and are canonical conjugate phase fields for j-th chain, and <f>± = (0 X ±(f> 2 )/y/2, 
v s is the spin wave velocity, and a is the lattice constant which we put 
to be 1. The mass m represents the spin gap, which originates from the inter-chain exchange 
interaction J±, i.e., 
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Antiferromagnetic J_l corresponds to the two-leg spin ladder, while the (strong) ferromag- 
netic J_|_ to the Haldane system. The magnetization density M(x), which is coupled to 
the magnetic field h, is given by the symmetric part as M(x) = d x (f> + (x) / \pK. It is known 
that both H + and H_ have the massive spectra when h — 0. Only H + is modified by the 
magnetic field h and becomes gapless, while the antisymmetric part remains unchanged. 
The coupling to the lattice distortion is given by 
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where Uj represents the dimerization of the j-th chain which modulates the intrachain ex- 
change interaction, while v the change of the interchain distance which modulates the in- 
terchain exchange interaction. We treat the fields Uj, v classically, which is justified be- 
cause the lattice is usually three-dimensional. It can be seen from eq.(2) that the mag- 
netic field will introduce solitons when the magnetic energy gain is larger than the soli- 
ton formation energy. To make it more explicit we transform H + and the coupling to 



the interchain modulation v into the Fermionic form by introducing the Fermion operator 
i>n, L (x) = (27ra o )- 1 / 2 exp[±^(0 + T 9+)] @: 
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This transformation is possible because H + describes the sine- Gordon model at f3 = An Hl8 
Equation (6) represents the Dirac Fermion with the mass m coupled to the lattice distortion 
v{x) with the Fermi energy controlled by the magentic field h. When h increases beyond 
h c \ = m, the Fermi level lies in the upper band of the Dirac Fermion. The magnetization M 
is represented by the Fermi wavenumber kp as kp = ttM where M is the normalized to be 
unity when all the spins are polarized. Then the generalized spin-Peierls transition occurs 
at the nesting wavevector Q = Ikp = 2ttM, i.e., 

v(x) = 2v cos(2vrMx + 5) (7) 

with 5 being some phase. The displacement is schematically shown in Fig. 2. First consider 
the case where the magnetization M is not too small and the Fermi energy ep is large 
enough compared with the transition temperature T c . A mean field treatment results in 
a second order phase transiton at T = T c = epe' 1 ^ N( - eF ^' y2 ^ with the specific heat jump 
AC/C(T C ) = 1.42 where C(T) is proportional to T. The validity of this weak coupling mean 
field theory can be examined as follows. The experiment has been done for h = 9.0T where 
the magnetization M = 0.3 p||9|,[IO|l . The gap m = O.SQmeV^ and the spin wave velocity v s is 
estimated from the intrachain exchange interaction J2 = 0.21meV as v s = |J2 — 0-33meV 
||. Therefore v s kp = irv s M = 0.3meV, which is about one-third of m. The Fermi energy ep 
measured from the bottom of the upper band is then Ep = 0.05meV = 0.5K. This estimate 
shows that the transition temperature T c is the same order as the Fermi energy ep, and the 
system is in the intermediate coupling regime. 

As the magnetic field is decreased to h c i the Fermi energy ep also decreases and the mean 



field treatment breaks down. Then we consider the zero temperature limit with changing 
the magnetic field. When the magnetization M is small and the off-diagonal matrix element 
7t>o introduced by the distortion is larger than the Fermi energy £p, many k points are 
connected by that matrix element. This is due to the parabolic dispersion which becomes 
important near the bottom of the band. As shown below the kinetic energy gain due to 
the distortion is given by AtiM^vq and vq is determined as vq = 2n^M. With this t>o, 
the matrix element is 2nj 2 M and equating this to the dispersion (v s k ) 2 /2m we obtain 
k = y/Anm^M/vs The number G of the reduced Brillouin zones within the region [—k , k ] 
is G = k /(nM) = ^4mj 2 /(7iv 2 M) and G » 1 in the limit of small M. In this case we 
can neglect the dispersion (v s k) 2 /2m in the diagonal matrix element compared with the 
off-diagonal ones, and the eigen-value problem can be approximately regarded as that for 
the tight binding wave with the wavenumber q discretized as q n = 2nn/G (n: integer). Then 
the eigen-value is given by 

e n = -2w 7 cos q n . (8) 

Therefore the lowest band (q — n — 0) with the energy — 2^v$ is occupied and the energy 
gain is —Ati^vqM. Then minimizing v 2 — Ati^vqM, we obtain v = 2n^M as mentioned 
above. The gap between the lowest ( n — 0) and the second lowest one (n — 1) is of the 
order of v^jG 2 oc M 2 and is much smaller than the matrix element t>o7(oc M) in the limit 
of small M. This happens when the magnetic field h is slightly above the critical field hfi 
without the lattice distortion and we define Ah = h — h^i ■ Based on the above discussion, 
let us consider Ah as a function of the magnetization M. Without the lattice distortion, 
Ah = (irv s M) 2 /(2m) for small M. With the coupling to the lattice, it changes into 

Ah = -2>yv + a(irv s M) 2 /m = -4ir>y 2 M + a(irv s M) 2 /m (9) 

where a is a constant of the order of unity. This is because the chemical potential at T = OK 
is at the middle of the gap between the lowest (n = 0) and the second lowest (n = 1) bands 
which is — 27t> + O ({v S M) 2 / m) . The relation (9) means that the magnetization M jumps 
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from to a finite value M = Mo ~ m^ 2 jv 2 with a first order phase transition. Summarizing 
the considerations above, the schematic phase diagram for spin ladder system coupled to 
the lattice is given in Fig. 1. The first order phase transition line terminates at a tricritical 
point to change into the second order transition line. 

There is no instability for the dimerization Uj, which corresponds to the standard spin- 
Peierls distortion, when the coupling constant g is small enough. This is because the operator 



sin v27T0j = sin-y/7r(0 + ± 0_) contains cos y^^- and sin ypUfy-. In ref. [|l8j it is discussed 
via the mapping the Ising model that the correlation functions for cos v^r^- and sin \pn<\>- 
decay exponentially. Then also the correlation function of sin \/2~K(j)j decays exponentially. 
This means that the energy gain due to the dimerization Uj is given by ~ g 2 u 2 for small g 
even at zero temperature, and the system is stable against the dimerization as long as g is 
small because of the elastic energy u 2 /2. Therefore, the lattice instability does not occur in 
the Haldane system where the interchain coupling J± represents the ferromagnetic Hund's 
coupling and is very strong (~ leV). The spin gap is due to the quantum fluctuation of 
the resulting 5 = 1 spins n, and the modulation of J± does not play any role. Then the 



magnetization process can be discussed by a model without the coupling to the lattice ||T9| , |20 
Another instability is of course the AF ordering as discussed in ref. The competition 
between the lattice instability and the AF is determined by the relative magnitude of the 
inter-ladder exchange interactions and the spin-lattice coupling, and the coexistence of these 
two is not expected for clean system |21 



We now discuss the experimental consequences of the above scenario, (i) The Bragg spot 
should be observed in the X-ray and/or neutron scattering at the wavevector Q = 2nM 
due to the modulation of the interchain distance, (ii) There appears a gap below T c and 
the specific heat should behave as C ~ e -^l T ; where 7^0 is the gap of the order of the 
transition temperature T c . (iii) The magnetization curve shows a hysteresis behavior at low 
temperatures. The experimental confirmation of these predictions is highly desired. 

In conclusion we have studied the lattice instability in the spin-ladder system when the 
spin gap is collapsed by the magnetic field. The lattice will distort to modulate the interchain 
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exchange interaction with the incommensurate wavenumber Q = 2irM. 
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Figure captions 

Fig. 1: Schematic phase diagram of the two-leg spin ladder system coupled to the lattice 
distortion in the plane of temperature T and the magnetic field h. The thick solid curve is the 
first order phase transition line while the broken one the second order transition between the 
undistorted and distorted states. The magnetization jumps across the first order transition 
line. 

Fig. 2: Schematic view of the lattice distortion which modulates the interchain exchange 
interaction. The wavenumber Q = 2nM changes as the magnetization M increases. 
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